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^ ' We develop a novel approach to gravity that we call 'matrix general rela- 

O ■ tivity' (MGR) or 'gravitational chromodynamics' (GCD or GQCD for quan- 

^ ! turn version). Gravity is described in this approach not by one Riemannian 

metric (i.e. a symmetric two-tensor field) but by a multiplet of such fields, 
or by a matrix-valued symmetric two-tensor field that satisfies certain condi- 
' tions. We define the matrix extensions of standard constructions of differen- 

tial geometry including connections and curvatures, and finally, an invariant 
functional of the new field that reduces to the standard Einstein action func- 
^ \ tional in the commutative (diagonal) case. Our main idea is the analogy with 

Yang-Mills theory (QCD and Standard Model). We call the new degrees of 
freedom of gravity associated with the matrix structure 'gravitational color' 
or simply 'gravicolor' and introduce a new gauge symmetry associated with 
this degree of freedom. As in the Standard Model there are two possibilities. 
First of all, it is possible that at high energies (say at Planckian scale) this 
symmetry is exact (symmetric phase), but at low energies it is badly broken, 
so that one tensor field remains massless (and gives general relativity) and 
the other ones become massive with the masses of Planckian scale. Second 
possibilty is that the additional degrees of freedom of gravitational field are 
confined within the Planckian scale. What one sees at large distances are 
singlets (invariants) of the new gauge symmetry. 
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1 Introduction 

Gravity is one of the most universal physical phenomena in Nature. At 
the same time it is also one of the most challenging problems in theoretical 
physics. In the Newtonian mechanics all gravitational phenomena are de- 
scribed by one scalar field (gravitational potential) subject to the Poisson 
equation. Einstein, when started to think about the nature of space and 
time, soon realized that it is impossible to develop a consistent relativistic 
theory of gravity with a scalar field and one needs a symmetric two-tensor 
field (which can be interpreted as the Riemannian metric of the space-time) 
subject to what is now called Einstein equations [Tj. Nowadays, Einstein 
General Relativity is accepted as a correct theory of gravitational phenom- 
ena at huge range of scales, from cosmological to the subatomic ones. In spite 
of the fact that General Relativity successfully describes all classical phenom- 
ena (with, maybe, few exceptions like singularities, dark matter, etc), so far 
it withstands all attempts to quantize it. In other words, we still do not have 
a consistent theory of quantum gravitational phenomena, that is phenomena 
at very small length scales (or high energies). It is expected that the gen- 
eral relativistic description of gravity, and, as the result, of the space-time, 
is inadequate at short distances. One needs new ideas to modify or to de- 
form General Relativity. There are many different proposals how to do this 
(string theory j2], noncommutative geometry j3j|, loop gravity |4j etc) but 
none proved to be the right one so far. An appealing idea is that all these 
approaches will be somehow related within one big unifying picture called 
M-theory 0. 

In this paper we start from the very beginning and carefully analyze 
the origin of the standard geometric interpretation of gravity. We show 
how this differential-geometric language can be generalized so that standard 
general relativity appears in a special commutative (or diagonal) limit. We 
propose that gravity should be described not by one two-tensor field but by 
a multiplet of tensor fields (or by a matrix valued two-tensor field) with the 
corresponding gauge symmetry incorporated in the model. Our approach 
should be contrasted with the "noncommutative extensions of gravity" on 
non-commutative spaces 13 E] ; it is also different from the model studied 
in 0. 
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2 Origin of Riemannian Geometry 

Let us recall the origin of Riemannian geometry. As a matter of fact its 
roots are in the theory of linear second-order partial differential equations of 
mathematical physics that describe wave propagation, in particular, light. 

Let M be a n-dimensional manifold without boundary. Our construction 
will be purely local, so it does not depend upon whether or not M is compact 
or noncompact. Let x denote points of M. We will be working in a small 
neighbourhood of a fixed point, say the origin, that can be covered by a single 
system of local coordinates x^, (/x = 0, 1, . . . , n — 1). 

Let / e C°°{M) be a real- valued smooth function on M (a source field) 
and let us consider an equation 

L^-f, (1) 

where Lp is an unknown function, and L : C°^{M) C°°{M) is a linear 
second-order partial differential operator with real smooth coefficients of the 
form 

L = a>'''{x)d^d, + lf{x)d^ + c{x) . (2) 

Most importantly, it is required that the matrix a'''^ e C°°{M) is a real 
smooth symmetric nondegenerate matrix, i.e. for any x & M 

o^^(x) = o'^^(x) , (3) 

deta'"'{x)y^O. (4) 

The determinant of the matrix a is a smooth function over M that does 
not change sign; it is either strictly positive of strictly negative. Since the 
matrix a is symmetric it has real eigenvalues. Let and n_ be the numbers 
of positive and negative eigenvalues of the matrix a (we assumed that there 
are no zero modes). Since a{x) is a smooth matrix valued function, the 
numbers n+ and n_ must be constant throughout the whole manifold, i.e. 
the signature of the matrix a is constant. This is because for an eigenvalue 
to change sign it must go through zero, but this would make the matrix 
degenerate. So, if the matrix is non-degenerate at every point and smooth, 
then the eigenvalues cannot change sign, and the number of positive and 
negative eigenvalues cannot change, i.e. they are stable. 

If n_ = 0, i.e. the matrix a is positive definite with the signature 
(-1- • — h), then the operator L is elliptic. Elliptic equations describe problems 
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in potential theory. All coordinates are physically of the same type, they are 
'space coordinates' and M is the physical space. 

If n_ = 1 then the matrix a has signature ( — !-■■■+) and the operator 
L is hyperbolic. Then there is one coordinate that is very different from the 
others. This is the coordinate along the negative eigenvector of the matrix 
a, and it is called 'time coordinate', t — x'^, versus 'space coordinates', x\ 
{i = — 1). Hyperbolic equations describe propagation of waves. In 

this case the manifold M is called the spacetime. One usually assumes that 
the spacetime M has the following topological structure M = M x S, where 
E is a manifold (compact or noncompact) without boundary (a time slice of 
M) . Since we are concerned only with local questions such topological issues 
will not play any role. 

Further, we note that under the smooth diffeomorphisms 

x'" = x'^ix) , 
the matrix a'^'^ transforms according to 

which is exactly the transformation law of the compoucuts of a contravariant 
two-tensor of type (2, 0), i.e. a section of the bundle TM ®TM . This can be 
used to transform a hyperbolic operator L to the following canonical form 

L = -d'l + a'^{t, x)didj + b%t, x)dt + b\t, x)di + c{t, x) . (6) 

Important information is provided by the characteristics of a hyperbolic 
equation: they define the wave fronts that serve to describe the connection 
between waves and particles (geometric optics). They are also needed to 
find the short wave asympotics of the solutions of hyperbolic equations. The 
characteristics are a family of the level curves S{x) — C oi the solution of 
the Hamilton- J acohi equation 

a>'^{x){d,S){d,S)^Q. (7) 

This equation is closely connected with the Hamiltonian system (x,^) with 
the Hamiltonian 

H{x,^)^a^^{x)U., (8) 
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which leads to the Hamiltonian equations 

2a^^^x)^, (9) 

-d,a'^^{x)Up. (10) 

The spacetime has a causal structure defined as follows. Let the initial 
conditions be x{0) = Xq and ^(0) = ^o- Then for a fixed Xq and varying the 
tangent lines to the trajectories with H{xq, ^q) = define a causal cone C{xq) 
at the point xq (usually called the light-cone). The causal cone separates 
M in two regions: the interior T{xo) of the cone that can be called the 
causal set (consisting of the trajectories with H(xo,Co) < 0. so called time- 
like trajectories) and the exterior S{xq) of the cone {causally disconnected 
set) (consisting of the trajectories with H{xq,^o) > called the space-like 
trajectories), so that 

M ^I(xo)UC(xo)u£{xo) . 

The points in the causal set are causally connected with the point Xq, i.e. 
they can be connected by time-like trajectories. The causal set is divided in 
two parts: the absolute past I~{xo), and the absolute future I'^{xo) of the 
point Xq, i.e. 

2:(xo) =X-(xo)UX+(xo). 

The exterior of the cone is not causally connected with Xq. 

The study of linear second-order partial difi:erential equations (both el- 
liptic and hyperbolic) simplifies significantly if one introduces the machinery 
of Riemannian geometry. The transformation law of the matrix a under 
diffeomorphisms enables us to identify it with a Riemannian metric 

g""" = 9,. = ian~' , (11) 

where {a'^'^) is the inverse matrix. 

Then, one defines the canonical Levi-Civita connection (Christoffel sym- 
bols) 

r'^a/J = -g^'^ {dpg^a + dag^p - d^Qafi) ■ (12) 

The Levi-Civita connection defines in a canonical way a connection on all 
bundles in the tensor algebra over the tangent TM and cotangent T*M 



dxf" 
~~dt 

dt 
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bundles, i.e. covariant derivatives of tensors of all types. In particular, 
the Levi-Civita connection is the unique compatible torsion-free (symmetric) 
connection. In other words, the Christoffel symbols are the unique solution 
of the equations 

= T^Pa (13) 

^a9t,v = d^g^,u - ^^^,agxu - ^^uag^^x = . (14) 

The Hamiltonian system (fTUjl is nothing but the equation of geodesies 
of the metric g. The geodesies lying on the surface H{xo,^o) = are the 
nuU-geodesics that form exactly the light cone at xq. 

Finally, one defines the Riemann curvature tensor 

R^uaP = doX^uf} — dp^^va + ^^XoX^ufl — ^^\l3^'^ua , (15) 

the Ricci tensor 

RvP = R'^UpLf) , (16) 

and the scalar curvature 

R = g^'R.u . (17) 

It is worth mentioning that the metric can be defined purely algebraically. 
Since L is a second-order partial differential operator, for a scalar function 
5* (which we view here as the operator of multiplication by the function S) 
the commutator 

[L, S] = LS- SL (18) 
is a first-order differential operator, and the double commutator 

[[L,S],S] = 2a^''id,S)id,S) (19) 

is just a smooth function. Thus the Hamilton- Jacobi equation (jT)) can be 
written in the form 

[mS],S] = 0. (20) 
This enables one to define the metric in terms of the commutators. 



3 General Relativity 

So far we were studying some matter field in a given background metric 
(a tensor field). The dynamics of the matter field can be described by the 
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matter action functional 

S„,^t{^,g) = j dxC{ip,g) , (21) 

where dx = dx^ A dx^ A ■ ■ ■ A is the Lebesgue measure on M, 5*) is 

the Lagrangian density and is a region of the spacetime with a spacehke 
boundary dVt. 

In general relativity one identifies the metric g with the gravitational field 
and studies the dynamics of the gravitational field itself by considering the 
Einstein-Hilbert action functional (with the cosmological constant) 

s,rM = Jdx^Agl (R J ^ ' ^^^^ 

n an 

where \g\ = |det(7^j,|, G is the Newtonian gravitational constant, A is the 
cosmological constant, x are the coordinates on the boundary dQ, g is the 
induced metric on the boundary and K is the extrinsic curvature of the 
boundary dQ. Again, the surface term is written just for completeness, it 
will not play much role in the following. That is why we do not provide any 
more details on the definition of the induced metric, the extrinsic curvature 
etc. This is all standard material that can be found in standard references 

HI 

The standard variation procedure leads then to the matter field equations 

j^S^^,iip,g) = 0, (23) 

and the Einstein equations 



where 



Rt.u - -g^^uR + f/^.A = 8nG T^, , (24) 



r^'^ = j=^S^^,{^,g) (25) 



is the energy-momentum tensor. The action of general relativity is con- 
structed in such a way that it is invariant under diffeomorphisms. This leads 
to the conservation of the energy-momentum in the usual way. 
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4 Origin of Matrix Riemannian Geometry 

General relativity is constructed by using the following fundamental objects 
and concepts: 



1 

ii 
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iv 
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Event (a real n-tuple • • • , x""^)), 

Spacetime (the set of all events), 

Topology of spacetime, 

Manifold structure of spacetime. 

Smooth differentiable structure of spacetime, 

Diff'eomorphism group invariance, 

Causal structure (global hyperbolicity) , 

Dimension of spacetime (in low-energy physics n — A), 

(Pseudo)-Riemannian metric with the signature ( — h • — h). 

Canonical connections on spin-tensor bundles over the spacetime. 



One can deform general relativity by changing the meaning of different 
aspects of this picture. The most radical is probably the approach of non- 
commutative geometry when one replaces the basic structures like events, 
spacetime etc. The most conservative approach is to change just the metric 
part without changing much of the above. Our approach is rather a conser- 
vative one since it just changes the least fundamental notions in the above. 

The analysis of the previous sections clearly shows that the basic notions 
of general relativity are based on the geometrical interpretation of the hy- 
perbolic wave equation that describe propagation of some fields without the 
internal structure, (in particular, light), that could transmit information in 
the spacetime. At the time of creation of general relativity the electromag- 
netic field was the only field that can be used for such purpose. This is still 
true for the macroscopic phenomena. However, at the microscopic distances 
this role of the electromagnetic field (photon) could be played by some other 
gauge fields (say, gluons and other vector bosons) that, together with the 
photon, form a multiplet of gauge fields with some internal structure. That 
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is why to repeat the Einstein analysis one has to consider a hnear wave equa- 
tion for such fields, i.e. instead of the scalar equation (j2)) we have a system 
of linear second-order hyperbolic (wave) partial differential equations. This 
would cardinally change the standard geometric interpretation of general rel- 
ativity. Exactly in the same way as a scalar equation defines Riemannian 
geometry, a system of wave equations will generate a more general picture, 
that we call Matrix Riemannian Geometry (MRG). 



4.1 Hyperbolic Systems 

To be precise, let F be a smooth A^- dimensional vector bundle over M, V* be 
its dual and End {V) ~ y x be the bundle of its smooth endomorphisms. 
Then the sections of the vector bundle V are represented locally by complex 
A^-dimensional contravariant vectors cp = (^Lp'^{x)), the dual vectors x ^ V* 
are complex covariant A^- vectors x = {xa{x)) and the sections X of the 
endomorphism bundle End {V) are represented by x complex matrices 
X={X^Bix)). 

We assume that the vector bundle V is equipped with a Hermitian fiber 
inner product that can be represented locally by 

{ij,ip) = i^EABv'', (26) 

where the bar denotes the complex conjugation and the matrix E = {E^b) 
defines the Hermitian metric, i.e. it satisfies the equation 

W=E, (27) 

where T denotes the matrix transposition. The Hermitian inner product 
provides a natural isomorphism between the bundles V and V* by 

{iP,^)= tiv{v®r), (28) 
where if)* G V* is the section dual to il) and tr y is the fiber trace. Locally 

rA = i^EBA=^I^ . (29) 
Similarly, we will also identify the bundles (End iV))* and End (V^) by 



(^,X^) = (X>,¥.), 



(30) 
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so that 

X* = E-^X^E. (31) 

Let (yj, / be smooth sections of the bundle V and let us consider the 
equation 

Lv = f (32) 

where L : C^iV) —>■ C°°{V) is a linear second-order partial differential 
operator of the form 

L = a^"'{x)d^d^ + 6^(x)a^ + c(x) (33) 

with endomorphism-valued smooth coefficients, i.e. 

^ (a/^-^^(^)) , = [lf^^{x)) , c = {c^b{x)) . (34) 

As in the scalar case a^'^ must be symmetric in the vector indices 

a^"" = a''^' . (35) 

Further, we assume that the operator L is formally self-adjoint with respect 
to the fiber inner product and some measure /i, which means that the com- 
ponents of the matrix a must be self-adjoint 

[a"")* = a^"" . (36) 

There are also some conditions on the lower order coefficients b and c but 
they will play no role in the subsequent discussion. 

Note that the matrix a transforms under diffeomorphisms as a contravari- 
ant matrix- valued two-tensor, (more precisely a section of the bundle TM ® 
TM End {V)), i.e. exactly as in eq. (0). However, in the matrix case one 
cannot, in general, put the operator L in a form like (jHI) by choosing the 
coordinates. We can also consider the gauge transformations 

V{x) — U{xMx) , (37) 

a'^^(x) — > U{x)a^''{x)U{x)'^ , (38) 

where U{x) is a smooth nondegenerate matrix-valued function. 
The leading symbol of this operator is 



aLmx,0 = -H{x,0 



(39) 
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where ^ G T^M is a cotangent vector and 

H{x,0 = a^''{x)U,. (40) 

Obviously, the matrix H{x,^) is self-adjoint 

H* = H (41) 

and, therefore, has real eigenvalues. We will assume that its eigenvalues 
hi{x,^), {i = 1, . . . , s), are distinct or have constant multiphcities di. 
It is clear that they are homogeneous functions of ^ of degree 2 

/i,(x,AO = A%(x,0. (42) 

Further, the eigenvalues are invariant under the gauge transformations ()38p 
and transform under the diffeomorphisms as 

h[ix',0 = h.ix,a, (43) 

where 



C = . (44) 



dx°' 

The operator L is elliptic if for any x E M and any ^ 7^ the matrix 
H{x,^) is nondegenerate, i.e. 



detvH{x,Oj^0, (45) 

so that all eigenvalues hi{x,^) are non-zero. 

The operator L is (strictly) hyperbolic at x in the direction u G T*M if 

det vH{x, p) = det v{a'"''{x)u^u^) (46) 

(i.e. the matrix a^^{x)v^Vy is non-degenerate), and for any cotangent vector 
^ 7^ not parallel to v all the roots of the characteristic equation 

det vH{x, e + Az/) = det y [X^a'^^'u^u, + 2\a^''^^v, + a^^'i^i,] = (47) 

are real (and distinct); clearly there are 2N roots. 

The same condition can be stated as follows: the operator L is (strictly) 
hyperbolic at x in the direction u eT*M if 

hi{x,u)j^O, (z = l,...,s) (48) 
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and for any cotangent vector 7^ not parallel to v each characteristic 
equation 

/ii(a;,^ + A//) = (49) 

has exactly two real distinct roots A^(a;,^) . 

The cotangent vector v[x) defines a smooth one-form on M and the roots 
of the characteristic equation are smooth functions on M. Let us denote 



V and choose tu^ . . . , {jJ^~^ G T*M so that c<j", (a = 0, 1, . . . , n — 1), is 
a basis in the cotangent bundle. Then the operator L is (strictly) hyperbolic 
if at any x G M none of the one- forms uj^, . . . jUj"^"^ defines a hyperbolic 
direction. 

The system of hyperbolic partial differential equation describes the prop- 
agation of a collection of waves. Similarly to the scalar case the operator L 
generates the causal structure on the manifold M as follows. First, we define 
the characteristics of the matrix hyperbolic operator L. The Hamilton- Jacobi 
equations and the Hamilton equations have the form 

det v[H{x, dS)] = det y [a^''{x){d^S){d^S)] = , (50) 
= 2tr,,-|A'(2:.e)a'-(a:)]5„ (51) 

= -tvv[K(x,i)dX''(=')](A (52) 

where 

K{x,0 = detv[Hix,0]H-\x,0- (53) 
Note that K{x,^) is polynomial in the matrix a^'^ . 

The Hamilton- Jacobi equation (j5U|) has then as many solutions as the 
number of eigenvalues. Each eigenvalue defines a Hamiltonian system of its 
own, i.e. a Hamilton- Jacobi equation 

hi{x,dS) = 0, (54) 

and Hamilton equations 

dx^^ d 



dt d^^ 
d^f, d 



dt dx^^ 



h,{x,0 (55) 
hiix,0 (56) 
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These equations define different trajectories for each eigenvalue. These 
trajectories can be identified with the geodesies of some Riemannian metrics. 
The trajectories with tangents on the surface det vH{xo, ^q) = are the null 
trajectories and define the causal cones. In general, they are different. So, 
instead of a unique cone one gets s causal cones Ci{xo) (recall that s is the 
number of different eigenvalues of H{x,^)). Each cone defines a causal set 
Xj(xo) consisting of the absolute past X~{xq) and the absolute future I^{xq) 
as well as the exterior of the cone Si{xo). Since the points in all causal sets 
are causally connected with the point xq, i.e. there is at least one time-like 
trajectory connecting those points with xq, we define the causal set as the 
union of all causal sets 

s 

I{xo) = [jliixo), (57) 

1=1 

similarly for the absolute past and the absolute future 

s 

2:±(a;o) = U^^K), (58) 

i=l 

The causally disconnected set is defined as the intersection of the exteriors of 
all causal cones 

s 

S{xo)=f]S^{xo). (59) 

i=l 

With these definition we have the standard causal decomposition 

M = I{xo) U dI{xo) U £{xo) 

^ 2:~(xo) ux+(xo) uaj(xo) U6:(xo) . (60) 

Since the causal cones vary from point to point, the structure of the causal 
set is different at different points. Such a picture can be interpreted as a 
"fuzzy light-cone." 

We see that in the matrix case the operator L does not define a unique 
Riemannian metric. Rather there is a matrix-valued symmetric 2-tensor field 
a'^^. We can decompose it according to 

where I is the identity endomorphism, k is a deformation parameter, 

^^^ = 4trya'^^ (62) 
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and (f)'^'^ is the trace-free part 

try0''" = O, (63) 
so that the matrix H{x, ^) becomes 

H{x,0 = ^9^''ix)^4'^ + K(P>^'{x)^,^,. (64) 
We can also introduce the whole family of 2/i;-tensors 

^Mii^i-M.i^fc = ^ tr va''"'' ■ ■ ■ a^"^"'' , (65) 

which contain the information about the matrix a'^'^. Similarly, the informa- 
tion about the eigenvalues of the matrix H{x, ^) is encoded in the traces 

2=1 

This can be easily evaluated as a power series in the deformation parameter 
^tTvH^{x,0 = {9'''Uuy (67) 

It is worth noting that, in general, the matrix g^^ is not necessarily in- 
vertible, even in the elliptic case. Although in the weak deformation limit, 
i.e. K — >■ 0, it must be nondegenerate, it is not necessarily so in the strongly 
deformed theory for large n. Thus, the matrix g^" which plays the role of 
Riemannian metric in the commutative limit loses this role in fully noncom- 
mutative theory and can be singular or even zero. 



5 Deforming General Relativity 

Gravity will be described by new dynamical variable a^'^ . Our final goal is 
to construct a diffeomorphism-invariant functional S{a) from the matrix a^^ 
and its first derivatives. Notice also that g^^ , in fact all tensors gt^^"'t^^'= ^ 
do not change under the gauge transformation ()38|) . Our main idea is to 
promote this symmetry to a universal local gauge symmetry. 
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That is why we need an invariant action functional of this field 

'S'matr-grav(a) = j dx C{a, dtt) , (68) 

where C{a, da) is the Lagrangian density. This functional should be invariant 
under both the gauge group as well as the group of diffeomorphisms. In 
infinitesimal form these transformations read 

S^a"" = [u, a""] , (69) 

5ga^- = ^^dxa'"' - a^'^dxC - a'^dxe , (70) 
where u is an element of the algebra of the gauge group 

= io (71) 

and ^ is the infinitesimal coordinate transformation 

6^x>^ = -e{x) . (72) 

The action for the matrix gravity should reduce to the standard Einstein- 
Hilbert functional (j22p in the commutative limit k (or when 0^^^ = 0). 

First of all, we need a measure, i.e. a density /i(a) that does not depend 
on the derivatives of a and transforms under diffeomorphisms like 

IJ,'{x') = fi{x)J{x) , dx' fi' {x') = dx fj,{x) , (73) 

where 

J{x) = det 



dx'^{x) 
dx'^ 



(74) 



or, in infinitesimal form, 

6^fi = 0, 5^fi = ed^fi + fido.r = dMn- (75) 

As a guiding principle we will require the correct commutative limit 

/^=V^ + 0(/€), (76) 

where l^^l = [a det g'^'^)~^ and a = 1 in the elliptic case and a = —1 in the 
hyperbolic case so that yU > 0. 



15 



We can construct a good candidate for the measure as follows. Let F{x, ^) 
be a scalar function constructed from the matrix a^^{x) and ^ that is invariant 
under gauge transformations and transforms under diffeomorphisms like 

F\x',0 = F{x,O, (77) 

where 



dx' 
^x'^' 

For example, we can define 



C = -^ie. ■ (78) 



F(x,0 = ^trv$(//(x,0) = ' (79) 

where $ is a positive function of single variable such that it decreases suffi- 
ciently fast as — >■ oo. Then the function 

^(x) = j diF{x,0 = j d^^tTy^{H{x,0) (80) 

transform as ()73|) and can serve as measure. The choice of the function $ 
should guarantee the convergence of this integral. Note that this choice is 
obviously not unique. 

We can write the measure in the form 

1 , , 

M=^trvP, (81) 

where p is a matrix-valued function that transforms like density, i.e. 

p\x') = p{x)J{x) . (82) 
For example, p can be defined by 

p{x) = / (i^$(x,0 • (83) 



Another good candidate for the measure can be obtained as follows. Let 

^ = -,e,,...,^e,,...u„a''''''---a^"''-, (84) 
nl 
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where e is the standard completely antisymmetric Levi-Civita symbol. By 
using the symmetry properties of the matrix a^'^, (|35|). (j36j) . one can prove 
that the matrix i/j is self-adjoint 

ij* = ijj. (85) 

We will require that this matrix is nondegenerate. Then the matrix ilj*ifj is 
positive definite so that we can define p by 

p = iri^)-'^'- (86) 

Obviously, the matrix p is self-adjoint, p* = p, and positive definite, p > 0, 
so that the measure is positive p > 0. 

In the commutative limit both of the above definitions lead to the same 
Riemannian measure (fTUI) . In the weak deformation limit we obtain 

p = v^|l + K'^try(0'^>"^)(ci(7^„^,^ + C2g^,ugap) + 0{k^)^ , (87) 

where g^^, = {g^'^)~^ is the inverse matrix and ci and C2 are some constants 
depending on the choice of the measure. 

Thus we obtain the simplest zero order term in the action 

A f 



n 

A 



j dx a/h|i + f^^-j^ trv{(l)'"'<p°''^){cig^,agui3 + C2gi,ugaf3) 
n 



8ivG 

+OiK')j, (88) 

where A is the 'cosmological constant' and G is the gravitational constant. 
This functional is obviouly invariant under both the gauge transformations 
and the diffeomorphisms. 

The dynamical functional Si{a) must depend on the first derivatives of 
a. Assuming that it is local, it must be quadratic in the first derivatives and 
non-polynomial in a, i.e. 

Siia) = -1- [ dxP^\^^,s^B^Dd^.a''^''Ad.a^'''c: (89) 
i67rG J 
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where P depends only on a but not on its derivatives. This can also be 
written as a sum of the terms like 

^'^""^ ^ 16^ / ^tiv J2P{a){da)Q{a){da) . (90) 

We see that our functional is reduced to nothing else but a generalized non- 
linear sigma model. In principle, one can find the form of the object P by 
studying the condition of diffeomorphism invariance in the perturbation the- 
ory in the deformation parameter k. Alternatively, which should be much 
easier, one can define tensors built from the matrix a and use them to con- 
struct the invariants. 

The total action in perturbation theory should have the form 

^matr-grav(o) = j dx^/\i\^{R - 2K) (91) 

1 



where RP'^'^^ is the curvature of the metric g, R is the scalar curvature, V de- 
notes the standard covariant derivative with the canonical symmetric connec- 
tion compatible with the metric g, and F, W and V are tensors constructed 
polynomially from the metric, g^i, and g"'^, and the Kronecker symbol. 

Therefore, in the weak deformation limit our model describes general 
relativity and a multiplet of massive tensor fields of spin 2 with mass param- 
eters of order A in unbroken phase. Depending on the potential terms there 
could be also massles fields as well as spontaneous symmetry breaking. This 
question requires further study. 



5.1 Matrix Geometry 

The Riemannian metric is nothing but the metric on the tangent bundle 
and general relativity is the dynamical theory of that metric. It can also be 
considered as an isomorphism between the tangent and cotangent bundles. 
Similarly, the quantity a'^'^ introduced in the previous section is an example 
of such an isomorphism of other bundles and our main idea is to develop a 
dynamical theory of such isomorphisms. 
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Let us consider the bundle Iso (T, T*) ~ Aut(T) of linear isomorphisms 
of the bundle T = TM ® V onto the bundle T* = T*M (g) V, where TM and 
T*M are the tangent and cotangent bundles and y is a vector bundle over M. 
The bundle V is not of the spin-tensor type since it is supposed to describe 
some internal structure of matter fields versus external one described by the 
spin-tensor bundles. The isomorphisms B : T —>■ T* can be identified with 
the sections of the bundle r*M®r*M®End {V) and and the isomorphisms 
A:T*^T with the section of the bundle TM O TM ® End (V^) . 

For a G TM ® TM ® End {V) to define an isomorphism, the equation 

a^V. = V (92) 

with any given i/j E T*M ® V , must have a unique solution 

= K^ij^" . (93) 

In other words, there must exist a unique solution, h G r*M®r*M®End {V), 
to the equations 

a'^"&.A = &MX' = (^^I. (94) 

This can be put in another form. Let Cj G T*M ®V be the basis in the space 
of one forms valued in V . Then the equation ()94|) has a unique solution if and 
only if the bilinear form Aij = (ej, Aci) is nondegenerate, i.e. det Aij ^ . 

Three remarks are in order here. First, assuming the nondegeneracy of 
the tensor g^^ one can obtain the solution of the equation (jMj) in form of a 
power series in the deformation parameter 

oo 
n=l 

which converges for small k. Second, even if the matrix g^'^ is degenerate (or 
even zero) the matrix b^i, might still be well defined, which would correspond 
to the limit k oo. Third, one can easily show that the matrix b^^, satisfies 
the equation 

b% = b.,, (96) 

but is not necessarily a self-adjoint matrix symmetric in tensor indices, more 
precisely, 

b% ^ b,, , b^, ^ b,, . (97) 
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In that sense the matrix a^'^ is nicer since it has additional properties (|35p. 

dnni)- 

That is why, one can use the matrices a'^'^ and to raise and lower 
tensor indices in the same way as the metric tensor. One must be careful 
however, since the matrix 6^1, is not symmetric and the matrices a'^'^ and 
do not commute for different indices. In particular, 

a'^ba^ifi'^ ^ ip" , b^aa^'ip. ^ ifia , (98) 

a^'^a'^'T^f, ^ a^^a^'^T^p , a^^^T'^ p ^ Kaa^^T'^p . (99) 
The equations ()94|) have important implications for the derivatives 

a'^^d^bpx = -d^a'^f'bpx , b^pdy^ = -d^pa^^ , (100) 

and, therefore, 

d^^b^p = -bayd^a^'bsp , a^a"^ = -a'^W^b^sa'^ . (101) 

Now we introduce matrix- valued coefficients A'^a/s that transform like the 
connection coefficients under the diffeomorphisms, i.e. 

, dx'f^ dx^ dx^ dx'^^ d'^x" 

Let be the tensor bundle of type {p, q) . We define a linear map 

V:TP(^V-^ TgP+i ® V (103) 

by 



i=l 

This map is well defined as an operator between the tensor bundles valued 
in V. 

Now it is not difficult to construct the matrix curvature and the matrix 
torsion. For a ip & T*M we compute 

{'D\)^,a - {V^v)u^.a = -Tl\,.u^x + T\,,{V^)^^ , (105) 
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where 

afxv d^-A. au di/^ -\- ^A. f}^^^ av «^ f5v-^^ a^i i (106) 

Next, we can define the matrix Ricci tensor 

= 7e%a. (108) 

and, by using the tensor a^^, even a matrix- valued scalar curvature 

7^ = a'^'n^^ . (109) 

Note though, that because of the noncommutativity the definition of the 
scalar is not unique . 

Now we need to fix the connection, i.e. to relate it somehow to the 
tensor a. To fix the connection we should impose an invariant compatibility 
condition. We will impose the compatibility condition in the form 

a^a"^ + ^"am«^^ + ^^Am«"^ = . (110) 
The solution of this equation is 



2 



where S is an arbitrary matrix valued tensor satisfying the symmetry condi- 
tion 

S^P- = -S'^"P . (112) 

In general, this connection is not symmetric. In the commutative limit it 
reduces to 

^"am = r^^I + ^ {S\^ + + 5a/) + 0{k) , (113) 

where F";^^ are the Christoffel coefficients of the metric g^^'^ and the indices of 
the tensor 5" are lowered with the metric g. Therefore, the torsion is reduced 
to the tensor S, i.e. 

r\^^s\^ + o{K). (114) 
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So, to have the torsion-free theory we have to set S = 0. However, in the 
noncommutative case the connection is not symmetric even if S" = 0. 
It turns out that the compatibihty condition of the form 

df^b^p - A\^bxp - A^p^b^x = , (115) 

which looks simpler on a first glance, is, in fact, more complicated. Because 
the matrix b^^ is not symmetric this equation cannot be solved in a closed 
form. For example, the solution of this equation does not have the familiar 
form of Christoffel coefficients 

A\f3 ^ ^ {d^bp^ + dpb^^ - a^^ . (116) 

The equation ()115|) can be solved only within the perturbation theory in the 
deformation parameter. We could also have used the symmetric part of the 
matrix b, i.e. b(f^,y), but the the inverse matrix would not be symmetric, so 
this does not simplify the solution after all. 



6 A Model of Matrix Gravity 

By using the matrix curvature we can now construct a simple generalization 
of the standard Einstein- Hilbert functional (with cosmological constant). We 
define 

5matr-grav(a) = j rfo; ^ tr y p (a^'^7^%„, - 2A) . (117) 

n 

This functional is obviously invariant under global gauge transformations 

a{x) Ua{x)U~^ . (118) 

One can easily make it local gauge symmetry by introducing a Yang-Mills 
field B valued in End (V) and replacing the partial derivatives in the defini- 
tion of the connection coefficients and the curvature by covariant derivatives 

d,-^d^ + [B„ ■ ] . (119) 

Thus we finally obtain an invariant action functional 

^matr-grav(a,i3) = J ^ tr yp | (a"'^7^%„, - 2A) 

n 

-^a'^^^M^a"''-^/?-} . (120) 
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where e is the Yang-Mills coupling constant, p is defined by eq. flHH|) or (|HH|1 



7^' 



8 +\E } - d -\E A^ } 



;i2i) 



-'Act 



^^apCT_^^pCTa_^CTapl, 



pu 



(122) 
(123) 



We should also impose an additional constraint on the tensor 5", for example, 
just put it to zero, 5" = 0, since it is not a dynamical degree of freedom. 

This functional describes the dynamics of the 'matrix metric field' a and 
the Yang-Mills field B. It is invariant under the diffeomorphisms 



V, 'T*^ — if^ ( ir\ 

and the local gauge transformations 

a{x) — * U{x)a{x)U~^ (x) 



U{x)B^{x)U-\x) - {d,,U{x))U-\x) 



(124) 



(125) 
(126) 



We could assume that the gauge group is a Lie group G, say a compact 
simple Lie group like SU{N), that B takes values in the adjoint representation 
of the Lie algebra of G and a takes values in the enveloping algebra generated 
by the Lie algebra. In the commutative limit this reduces to the standard 
Yang-Mills fields coupled to gravity. 

Now we need to introduce interaction of the fields a and B with matter 
fields in such a way that will lead to the spontaneous breakdown of the 
gauge symmetry, so that in the broken phase in the vacuum there is just one 
tensor field, which is identified with the metric of the space-time. All other 
tensor fields must have zero vacuum expectation values. In the unbroken 
phase there will not be a metric at all in the usual sense since there is no 
preferred tensor field with non-zero vacuum expectation value. Alternatively, 
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one could expect the gauge (grav-color) degrees of freedom to be confined 
within the Planck scales, so that only the invariants (grav-white states) are 
visible at large distances. For example, at large distances one could only 
see the diagonal part g^'^ = tr ya'^'^, which defines the metric of the space- 
time at large distances and the gauge invariants like tr vT^'^ap-y , etc., which 
determine in some sense the curvature of the spacetime. 

One should stress that this model is a realization of a consistent interac- 
tion of tensor fields with gravity, which usually constitutes a problem. Notice 
that this model is, in fact, nothing but a generalized sigma model. So, the 
problems in quantization of this model are the same as in the quantization 
of the sigma model. 

7 Conclusion 

A careful analysis shows that the origin of Riemannian geometry in general 
relativity lies in the theory of the wave equations. We propose to replace 
a single wave equation by a system of wave equations at small distances. 
This brings completely new geometrical picture in the theory of gravitational 
phenomena. Instead of one Riemannian metric we have now a matrix- valued 
tensor field a^'^, which is the main dynamical field describing gravity. We 
also introduced a new gauge symmetry which is responsible for mixing the 
new degrees of freedom and a new Yang- Mills field B^. We constructed a 
second order action functional that describes the dynamics of the fields a and 
B and is invariant under the diffeomorphism and the new gauge transforma- 
tions. This functional may be viewed as a "noncommutative deformation" of 
Einstein gravity coupled to a Yang-Mills model. We introduce a deformation 
parameter k such that the theory has a "commutative limit" k — > 0. In 
the weak deformation hmit our model describes Einstein gravity, Yang-Mills 
fields, and a multiplet of self-interacting two-tensor fields that interact also 
with gravity and the Yang-Mills fields. We speculate that the new degrees 
of freedom could only be visible at Planckian scales, so that they do not 
exhibit themselves in the low-energy physics. However, the behavior of our 
model at higher energies should be radically different from the Einstein grav- 
ity since there is no preferred metric in the unbroken phase, when the new 
gauge symmetry is intact. It would be very interesting problem to study 
simple solutions of this model, say a static spherically symmetric solution, 
which would describe a "non-commutative black hole" . There are reasons to 
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believe that this model could be free from singularities. 

Of course, before one can take the matrix gravity seriously various im- 
portant (and interesting) questions have to be clarified, in particular: i) 
classical (commutative) limit, ii) quantization. Hi) semiclassical approxima- 
tion, iv) renormalization, v) spontaneous symmetry breaking, and, finally, 
vi) Planck confinement. 

We would like to make a couple of final remarks. First of all, to study this 
model in the one-loop approximation would require new methods since the 
partial differential operators involved will not be of the so-called Laplace type, 
i.e. they will not have the scalar leading symbol. Most of the calculations in 
quantum field theory were restricted so far to the Laplace type operators for 
which nice theory of heat kernel asymptotics is available. However, the study 
of heat kernel asymptotics for non-Laplace type operators is quite new and 
the methodology is still underdeveloped. For example, even the first heat 
kernel coefficients [Aq, Ai and A2) needed for the renormalization in four 
dimensions are not known in general. For some progress in this area (the 
calculation of ^1) see jTUl ITT] . 

Instead of the second-order operator L we could start from a first-order 
operator. Let D be a first order formally self-adjoint partial differential 
operator 

D = a^{x)d^, + I3{x) , (127) 
where a is an anti-self-adjoint 

smooth matrix-valued vector field and /? is a matrix valued scalar field, acting 
on the smooth sections of the vector bundle V . The square of this operator 
defines a second-order operator 

L = D^ = a^"'{x)d^d^ + b^{x)d^ + c{x) (128) 

with 

a^'a'' + a'^a" = 2a'"' . (129) 

We stress here once again that we do not assume that a'^'^ = or even 
a^^ = g^'^E with some automorphism E. This enables one to repeat the 
whole construction in such a way that the dynamical variables of gravity will 
be the generalized 'Dirac matrices' instead of a'^'^ . 
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Another possibility is, following the approach of [HIEIIZI; to extend our 
model to the spaces where the coordinates do not commute, which is achieved 
by replacing the usual products by the (noncummutative) star (Moyal) prod- 
uct 



/(x.)*9W=exp(l.-^^ 



/(x + u)g{x + v) 



u=v=0 



(130) 



where 9 is an antisymmetric tensor. 
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